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Abstract 

The CMB bispectrum generated by second-order effects at recombination can be calculated ana- 
lytically when one of the three modes has a wavelength much longer than the other two and is out- 
side the horizon at recombination. This was pointed out in [1] and here we correct their results. 
We derive a simple formula for the bispectrum, = — (1/6+cos 20)-(l— l/2-dhi{l 2 s Cs)/d\nls)-, 
where Cs is the short scale spectrum and 9 the relative orientation between the long and the 
short modes. This formula is exact and takes into account all effects at recombination, includ- 
ing recombination-lensing, but neglects all late-time effects such as ISW-lensing. The induced 
bispectrum in the squeezed limit is small and will negligibly contaminate the Planck search for 
a local primordial signal: this will be biased only by /^l ~ — 0.4. The above analytic formula 
includes the primordial non-Gaussianity of any single-field model. It also represents a consis- 
tency check for second-order Boltzmann codes: we find substantial agreement with the current 
version of the CMBquick code. 



1 Introduction and main results 

The Planck satellite [2] is currently taking data and has as one of its main objectives the 
exploration of non-Gaussianities in the /nl ~ few range. The primary interest is of course in 
primordial non-Gaussianities generated in the early Universe, but for such small values of /nl 
all second-order effects connecting the initial conditions to the observed CMB are potentially 
relevant. The theoretical control of all these effects is still not completely satisfactory. Though 
the complete second-order equations have been studied, see for example [3J, SI 0, El [TJ, El [9J, [TU1 
HU [12l [131 EH El El H7j , getting the full bispectrum induced by second-order effects requires 
to resort to a numerical code. To our knowledge, the only publicly available code including 
all second-order effects is CMBquick [18J. This is hardly satisfactory as the complexity of the 
bispectrum calculation and the importance of the result would require a cross-check among 
different programs. 

However, as pointed out in [T] the full bispectrum calculation becomes quite easy in the limit 
in which one of the modes is out of the Hubble radius at recombination, while the other two 
have a much shorter wavelength. The reason for the simplification is intuitive: a mode which 
is out of the Hubble radius cannot affect any physical process, so that the effect of the long 
mode will only intervene in the way the short-scale 2-point function is eventually observed. The 
purpose of this paper is to review and correct the results of [I] and to compare the obtained 
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bispectrum with the full calculation of the second-order Boltzmann code CMBquick [18], the 
publicly available Boltzmann code created and maintained by one of the authors of this paper. 

To calculate the bispectrum we have to evaluate the temperature anisotropy at second order 
(the explicit calculation in the fluid approximation is reviewed in appendix [A]), in the presence of 
a short and a long mode, much longer than the Hubble radius at recombination. The long mode 
has no physical effect at recombination: it just acts as an unobservable coordinate redefinition, 
as we will explain in Section [2j while recombination takes place perturbed — at linear order — by 
the short mode. If the long mode were even longer than the present Hubble radius this would 
be the end of the story, as its presence would not leave any effect whatsoever. But for the 
calculation of the bispectrum, we are interested in the situation in which the mode eventually 
comes back into the Hubble radius and we have to compare different large portions of the sky 
with the same physics at recombination, but modulated by a long mode. 

The long mode gives rise to various effects that will be discussed in Section [3j First of 
all, it gives a long- wavelength modulation of the temperature so that, even if the short-scale 
perturbations have the same dynamics in the various regions, the perturbations are defined 
with respect to a different average temperature in different directions of the sky. This gives the 
same effect as a local primordial non-Gaussianity, with /j^ = —1/6, and it was missed in PQ{^] 
The long mode also gives rise to geometrical effects: the relation between a physical scale at 
recombination and a coordinate length will be modulated by the long mode and this will induce 
an angular stretching of the 2-point function in different regions of the sky. Thus, the effect 
of the long mode will be proportional to the derivative of the 2-point function with respect to 
the angular separation in the sky. Finally, lensing close to the last scattering surface will also 
displace the observed 2-point function and part of this effect will also depend on the relative 
orientation between the long and the short modes. 

We will get the final bispectrum in Section [4j We can write it as a function of the power 
spectra of the long and short modes, C\ L and Cj g , and the relative angle between them, 9, in 
the form 



B lLlBl8 = C lL C,J\. + (ivo*2»)[2 dln Vs C is 
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This formula is exact in the squeezed limit as it includes all second-order effects at last scattering 
(including the early ISW effect and the effect of anisotropic stresses). Additional contributions 
to the bispectrum will come at much lower redshifts and are therefore physically distinct and 
can, in principle, be observationally separated. When the Universe starts accelerating, a large 
ISW-lensing signal is induced pjJJ EOJ Ell EI] , while when structures begin to form a host of 
non-linear effects takes place, especially on short scales, which are far from Gaussian (see for 
example [231 El] ) • Notice that the formula above takes into account the correlation of lensing 
(induced by the long mode) with the temperature fluctuation at last scattering: this is important 
as there is a relevant cancellation [lj between the isotropic part of lensing and the redefinition 
of spatial scales discussed above. 

The comparison of our analytical result with the current version of CMBquick |18] is success- 
ful (Section [5]), with some residual discrepancy most likely coming from the numerics, especially 
in the treatment of neutrinos. Notice that for the code the check is far from trivial: though in 
the squeezed limit the understanding of the CMB bispectrum is much easier, in the code the 
result is not dominated by few simple effects. On the contrary, all second-order effects must 
conspire to reproduce the simple formula above. 

We will show that eq. ([I]) already takes into account all the (usually small) primordial non- 
Gaussianity of single- field models. From this point of view, it can be seen as a clear observational 



1 The presence of this term was shown in the limit where all the modes are out of the Hubble radius in 
[3] and was later understood as a consequence of the local redefinition of the average temperature in [3] . 
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prediction of all these scenarios and as an observational extension of the Maldacena consistency 
relation [251 HE W\ ■ 

The bispectrum that we obtain is rather small and unlikely to be a worrisome source of 
contamination for the forthcoming Planck analysis. This is quantified in Section [6j where we 
find that second-order effects at recombination would appear as a small bias to a search of a 
local primordial non-Gaussianity by /^l ~ —0.4. This implies that, in the squeezed region, 
where the primordial local signal is peaked, we have an analytical control of the background 
induced by second-order effects. Conclusions are drawn in Section [7} 



2 A long mode as a coordinate transformation 

In this section we will discuss how second-order perturbations, in the limit in which one of the 
modes is out of the Hubble radius, can be obtained simply by a coordinate redefinition of the 
first-order results. We will start showing that it is possible to rewrite locally a perturbed FRW 
metric as an unperturbed one by reabsorbing a very long wavelength mode — a mode outside the 
Hubble radius — with a coordinate transformation. 

We start from a linearly perturbed metric in Newtonian gauge and we consider only scalar 
perturbations. The metric reads 

ds 2 = a 2 (rj) [-(1 + 2<S> L )dr] 2 + (1 - 2$ rfSijdx* : dx j ] , (2) 

where &l and are very long wavelength modes, thus practically independent of the position. 
We consider the following coordinate transformation [28] 

fj = r) + e(ri), 

x { = x\l - A) , U 

where e is an arbitrary function of time and A is an arbitrary constant. This is the most 
generic coordinate transformation of an unperturbed FRW solution which leaves the metric in 
Newtonian gauge. In the new coordinates the metric takes the form ^ with potentials given 
by 

$L = 3>L-e ~Ue , (4) 
^ L = m L -\ + Ue . (5) 

Thus, we can reabsorb the long-wavelength mode by choosing e and A such that §l = and 
§ L = 0, i.e. 

q> L = e > + Ue, m L = \-%e. (6) 

As explained in [28] , not all the choices of e and A give rise to the long- wavelength limit of a 
physical mode in Newtonian gauge. Indeed, e and A must satisfy certain conditions imposed by 
those Einstein's equations that vanish in the limit ki —> 0, where ki is the wavevector associated 
to &l and ^l- These conditions are 

{%' -U 2 )v = + , (7) 

where v is the velocity potential, and 

$ = ^ - 8irG5a , (8) 

where 5a is the anisotropic stress, defined by writing the spatial part of the stress-energy tensor 
as Tij = gijp + didjda. Using eqs. §6§ and the first condition ([7]) one finds that v = — e and that 
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A simply coincides (up to a sign) with the comoving curvature perturbation £ = — ^ + Hv (q 
in the limit ki — > 0, 

A = -Cl- (9) 
Using the above equation the second condition, eq. Q, reads 

d + 2Ue = -Cl - 8nG5a L , (10) 

which can always be satisfied by a suitable choice of e. In particular, integrating this equation 
yields 

e = J a 2 (( L + 8irG5a L )dT] . (11) 

For instance, integrating the above equation in the absence of anisotropic stress, one finds 
for radiation dominance (a oc rj) e = — tj^l/S, and for matter dominance (a oc rj 2 ) e = —t]Cl/5. 
More generally, one can describe the transition from radiation to matter dominance using a{rf) = 
a eq (2ar] + a 2 rj 2 ), where a is a constant which fixes the normalization of rj, and a oq is the value 



of the scale factor at radiation/matter equality. In this case eq. (11) gives 



f( \r f(\— 20 + 15arj + 3a V 

e = vma , m = — 15{2+ar]) 2 ■ ( 12 ) 

Now we want to use the coordinate change above to study the second-order evolution of a 
short- wavelength mode k$ under the modulation of the long mode k^. In the squeezed limit, 
kg 3> ki, the physical effect of the very long wavelength mode simply amounts to a modifi- 
cation of the time and spatial coordinates of the short-wavelength mode |25| I26j. We start from 
a perturbed metric in the coordinates fj,x l , where the long mode has been reabsorbed in the 
background quantities, 



dsf = a\fi) -{1 + 2<f> s )dfj z + (1-2V s^ijdx'dZ 1 . (13) 

The potentials $5 and can be written in terms of the coordinates (77, x l ) using the coordinate 
transformation ([3]). They become (for this transformation in matter dominance see [29]) 

* s (fj, = *S + e^--\x^, 9 8 (ij, 5») = *s + e-^ - Xx^ • (14) 



By performing the coordinate transformation ^ on the entire metric (13), this can be written 
as a second-order perturbed one, 

ds 2 = a 2 (i]) [-e^drf + e - 2 %cfeW] , (15) 

with 

$ = $ s + $ x + e—^- + Clx^ , (16) 
or/ ox 1 

* = *s + y L + e^- + Clx 1 -^ . (17) 

Here we have dropped the tildes to simplify the notation. All the quantities labelled by a 
subindex S are to be computed at first order only. Note that there are only scalar modes in 
the metric. Indeed, vector and tensor modes generated by second-order scalar perturbations are 



2 In [28] and other references this variable is more commonly called 1Z. Here we use the notation of 
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suppressed in the limit of — > 0. Note also that we have used the exponential form for the 
metric (15) to reabsorb cross-coupling terms between long and short modes in a nice form. 

The coordinate transformation affects also other scalar quantities such as the temperature 
T, which transforms as 

rp rp . 9T S idT S . s 

T = T S + e— — + ( L x — -j . (18) 
or] ox 1 

Defining the temperature perturbation O as 

Tfas^Tfa) (1 + 6(17,3;*)) , (19) 

one obtains, using T oc 1/qJj 

9 = e s - He(l + Q S ) + e 9 ^ + ( L x^ . (20) 

For later purposes, it is convenient to give here also the effect of the coordinate redefinition 
on the fluid velocity. As the fluid four-velocity is a vector, under a coordinate change its 
components transform as 

u»(r ] ,x*) = —u»(f},x i ). (21) 
From the expression above we obtain, for the three- velocity of the fluid v i (u l = e*t//a), 

^ +£ M +C ^M (22) 

In section [5] we have checked numerically that these expressions accurately reproduce the second- 
order evolution in the squeezed limit. 



3 The effect of a long mode on the CMB anisotropics 



In the limit of instantaneous recombination and considering only matter dominance, the observed 
CMB anisotropics o bs are given, in the squeezed limit, by 



e ob s(ra) = [0 + $ + n • v + $0 + $ 2 /2 + (0 + *)n • v] (77 



veci -Tec ) 1 



(23) 



where n is the direction of the line of sight. On the right-hand side, is the intrinsic temperature 
fluctuations, <3? is the gravitational potential and v is the velocity of the baryon-photon fluid; 
these quantities are evaluated at the physical recombination time r] rec and position x rec . A 
detailed derivation of this formula can be found in appendix [A] (see eq. ( 74 ) ) . 

In order to study the effect of a long-wavelength mode on the CMB anisotropics, we will 
apply the change of coordinates ^ on each of the terms on the right-hand side of this expression. 
Despite the assumptions in the derivation of eq. (23), i.e. matter dominance and instantaneous 



recombination, the results that we will obtain are much more general. We will come back on 
this point at the end of this section. 



In matter dominance eqs. (|16|), (|20|) and (|22|) become 

5 



1 d$s 
3 in r] 



3 dx l 



= 



s 



4 + 



|Mi + e s ) + i^|^ 

3 3 a m 17 



5 jdQs 
3 dx l 



1 



— ^ 
3 d In r\ 



3 dx l 



(24) 
(25) 

(26) 



3 Note that the local term ~He®s would have been absent, had we defined the temperature pertur- 
bation in exponential form, i.e. T oc Te , such as for the potential 4> and 
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where we have used that in matter dominance Cl = —5&l/3. We can now evaluate the effect of 
a long- wavelength mode on o b s by applying these equations to replace and v on the right- 



hand side of eq. ( 23 ) . The right-hand side rearranges nicely and the temperature anisotropy can 



be cast in terms of observable quantities, as 



@obs(ra) = ©obs,s(n) + B obs ,L(n) + o bs,L(n) 1 + 
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5n ■ Vn) @obs,s(n) 



with 



<3>obB,s(n) = [&s + $s + n ■ v s ] (rj iec , x TCC ) , 



(27) 

(28) 
(29) 



(29). 



In eq. (27), the dependence of o bs,s(^) and o bs,L(^) on r/ mc can be read from eqs. (28) and 



In these expressions rj rec corresponds to the physical time of photon emission. At zeroth 
order r] Tec is just given by the unperturbed recombination time 77* and it is therefore space- 
independent. In the presence of a long-wavelength mode, recombination will not take place 
at the same coordinate time but at constant physical temperature. In matter dominance the 
coordinate transformation eq. ([3| becomes 



f} = rj (l + -$ L ) 



x l (l - 



Recombination takes place at constant fj = 77*, which then gives 

^?rec " 



1-^1* 



(30) 



(31) 



Similarly, also the physical time of observation r] \> s , which at zeroth order coincides with the 
unperturbed observation time 770, will be changed by the presence of a long- wavelength mode, 
Vobs = (1 - $l(z)/3)77o. 

The presence of the long-wavelength mode affects also the relation between the direction 
of observation h and the physical position at recombination x rec . In the absence of 3>£, x rec is 
given by the zeroth-order geodesic equation, x rec = h(rjo — 7?*), while in the presence of <3?l it is 
given by 



n 



1 



1 



'/<> 



+2n / <$> L {x)dr]-2 \ (r7-77*)V_L$ L (x)dr? , 

(32) 

we have rewritten 



as a solution of the first-order geodesic equation. Note that in eq. (|27 
^rcc • V^ rcc by using the zeroth-order geodesic equation, x rec = n(ryo — 77* 



Equation (27) satisfies an important consistency check [4j: a constant mode which is still far 
out of the Hubble radius today does not affect the observed temperature anisotropy. Let us see 
explicitly that when replacing $>l = C, where C is a time and position independent constant, 
this constant drops out from the physical observables. To do so, one has to realize that there 



are hidden second-order contributions in the first two terms on the right-hand side of eq. (27) 
First of all, a constant mode also redefines the average temperature observed in the sky, i.e. 



<T(n)> = (f (1 + 8 obS)i (n))) = f 1 + 



C 
3" 



(33) 



up to first order in C. The observed anisotropy must be defined in terms of this average 



temperature, see eq. (59), while equation (27) above is defined in terms of the unperturbed 
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temperature T. The variation of the average temperature changes the amplitude of small scale 
perturbations, 

ew(n) • = e obSi5 (n) ■ (i - • ( 34 ) 



This C-dependent contribution exactly cancels with the first term in parentheses of eq. ( |27| ), 
using o bs,L = C/3. For simplicity we did not include this effect of redefinition of the average 
temperature in eq. ( |27[ ); indeed, this is only relevant when a mode is still out of the Hubble radius 
today, while modes of interest for the CMB bispectrum do not change the average temperature. 
As discussed above, the presence of a long-wavelength mode perturbs the coordinate time of 



recombination, eq. (31). This introduces a second-order dependence on C through eq. (28) 



C d 

0obs,s(77rec,ra) = G obSi5 (r?* , ti) - — ■ g @obs,s(V*A) , (35) 



at first order in C. This cancels with the second term in parentheses in eq. (27), again using 

@obs,L = C/3. 



The first-order relation between h and the position of the photon emission x rec , eq. (32), 
introduces another second-order contribution. For a constant mode, the lensing term in eq. (32) 
vanishes, while the first two terms yield 

x T ec(n) = n (l - ^ j (770 -i]*)+n 2C (r/ - 77*) = n (l + j (r/ - 77*) . (36) 



This C-dependent term cancels with the last term in parentheses in eq. (27). The constant 
mode redefines the coordinates; however, the relation between the direction of observation and 
the physical point of emission remains unchanged. 

It is important to notice that the long-wavelength relation 

e ob s,L = ^l = -Ul, (37) 
3 5 

where 3>£, is the Newtonian potential during matter dominance, does not assume that recombi- 
nation takes place in exact matter dominance. It holds also taking into account the transition 
from a radiation dominated Universe and including the effect of anisotropic stress. This is easy 
to understand as the standard calculations leading to the Sachs- Wolfe formula do not need to 
be done at recombination, but still hold later on when the Universe reaches exact matter domi- 
nance (assuming that the mode is still out of the Hubble radius). This can be checked explicitly 
as follows. The intrinsic temperature fluctuation for a mode which is out of the Hubble radius 
can be obtained by transforming to the unperturbed coordinate 77. Given that the temperature 
scales as T oc l/o, it is given, independently of the background equation of state, by 

G = -Ue . (38) 

Taking into account the gravitational redshift and the (early) ISW effect, we have 

Gobs = -He + $ rcc + / dr] ($' + = -Tie - * rec + 2$ MD = C + 2$md = -^C, (39) 

where we have used that $md = ^MDj eqs. Q and Q and the matter-dominance relation 



£ = — 5$/3. This implies that, though our derivation of eq. (27) holds in matter dominance, the 



first and third terms in parentheses — which are the ones that enter in the final bispectrum — will 



remain the same also when we depart from exact matter dominance. Indeed, eq. (37) implies 



that the relationship between the long-wavelength temperature fluctuation and the spatial re- 



definition of the coordinates, i.e. remains the same, and this fixes the factor —5 in (27). 
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Moreover, the factor of 1 in parentheses remains the same. It just depends on the change in the 
average temperature given by the long mode and it is required by the validity of the previous 
argument about the cancellation of an homogeneous mode. 

Another approximation that we have made so far, which can be relaxed without affecting the 
final result, is that temperature perturbations for the short modes are created instantaneously 
at recombination. Of course recombination is not instantaneous and moreover perturbations 
are also generated by the early ISW effect. However, it is easy to see that our results are not 
changed, because the long mode is taken to be out of the Hubble radius during the whole period 
when recombination happens and matter dominance is finally reached. For example, the fact 
that recombination is not instantaneous is encoded in the presence of a window function g which 
weights the source term S, schematically 

d V g(v)S(v). (40) 

Now the long mode will modify both the source term and the window function, through a time 
redefinition: the perturbed source term will be equal to the unperturbed one once written in 
the tilded coordinates S(r]) = So(fj) and the same will hold for the window function (which 
has a different transformation property, as it is a probability density in r/): g(r/)dr] = go(fj)dfj. 
The integral remains the same and this parallels the cancellation discussed after eq. ( |35| ) in the 
instantaneous limit. The only effect of the time redefinition is the geometrical effect of equation 



(32). Along the same lines, one can show that also the Silk damping effect is taken into account. 



4 The bispectrum in the squeezed limit 

In this section we derive the bispectrum in the squeezed limit, eq. 0, starting from equations 
(27) and (32). The idea is to calculate how a long mode affects the short-scale 2-point function 
and from this to get to the bispectrum, similarly to what done in the derivation of the consistency 
relation for the primordial 3-point function |25} l26t 127] . For simplicity we work in the flat-sky 
approximation where the angular direction can be written as h = (1,777,) and the multipole 
expressions are simply Fourier transforms with respect to the m's. 



Let us examine eq. ( 32 ) . The last term gives lensing and we are going to discuss it separately 
below. The first term in brackets does not contribute to the 3-point function as it depends only 
on the Newtonian potential at observation. The second term in brackets changes the distance to 
the last scattering surface, but its effect is suppressed by 77*/t7o and can thus be neglected. The 



second term in eq. (32), which involves the integral of $ along the line of sight, is the so-called 
Shapiro time delay [30]. This is suppressed because the integral tends to average out along 
the line of sight, unless the mode is of the order of the present Hubble scale: this effect can 



be neglected if we disregard the very first multipoles. In eq. (27) the time derivative combines 



with the change in the recombination time, eq. (31), as we discussed in the previous section: 
one is left with the geometrical effect that light is emitted at a different coordinate time, but 
this is again suppressed by t/^/t/o- The same geometrical suppression affects the radial part of 
the spatial derivative. Therefore, the effect of the long mode on the short-scale 2-point function 
only amounts to a stretching perpendicular to the line of sight and to the constant term, the 
first in parentheses in eq. (|27l). This gives 



(@obs("i2)Oobs(rn<3)>L = (G obs (m 2 )e o bs(m3))o+G obS;Z/ (mi) [2 - 5m • V m J (G obs (m2)G obs (m 3 )) , 

(41) 

where the long-wavelength mode is evaluated at the midpoint rh\ = (777,2 + mz)/2. Notice that 
the short scale 2-point function only depends on the distance rh = 777,3 — m 2- The factor of 2 in 
the parenthesis takes into account that the long mode affects both G obs (777,2) and G ob s(w-3). 
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To get to the 3-point function one has to correlate the expression above with the long- 
wavelength temperature and then take the Fourier transform. These steps are the same as in 
the derivation of the 3-point function consistency relation [251 ESI EZ] and one gets 

Note that the long-mode power spectrum C\ L should be understood only to contain the Sachs- 
Wolfe effect and not the late ISW. The same disclaimer applies to the analogous equations 
below. 

To this we have to add the contribution due to lensing, which comes from the last term of 



eq. (32). In this paper we are not dealing with late-time effects, but only with effects which 
occur close to the surface of last scattering. This means that we only have the contribution 
coming from the correlation of the lensing potential due to the long mode with the temperature 
anisotropy created at last scattering, and not the (much bigger) effect due to the correlation 
between lensing and (late) ISW effect |19|, [2Q| [2~Tl [22] . In the limit in which the lensing mode 
is out of the Hubble radius at recombination, the lensing bispectrum is easy to calculate and is 
given by [I] 



B l L lsh ~ GC Il C Is 



2 cos 26* - (l + cos26>) 



d\n(l 2 S C ls ) 



(43) 



dials 

where 9 is the angle between the long and short modes. It is important to stress that, although 
this result was obtained assuming perfect matter dominance and therefore neglecting the role of 
the early ISW and the anisotropic stress, the final result does not change when these effects are 



included, similarly to what happens with eq. (37). Indeed the correlation between temperature 



and lensing potential peaks at the time when the long mode comes back into the Hubble radius 
(before that it is suppressed by the lensing window, later it loses correlation with the temper- 



ature fluctuation [4]). This implies that (43) is not modified by recombination taking place 
when radiation is not completely negligible, because sufficiently long modes come back into the 
horizon in exact matter dominance and the relationship between the Newtonian potential and 



the temperature is unchanged, see eq. (37). 



The two equations (42) and (43) have a similar structure and can be put together to give 
B lLlsls = C lL C ls (l + 6 cos 29) ( 2 -^f^ 1 ) , (44) 

which is the main theoretical conclusion of the paper. In the squeezed limit, a local primordial 
non-Gaussianity would give a bispectrum of the form — 12f^Ci L C[ s . Therefore, this bispectrum 
corresponds to 

' (1 + aco.*)(i-i!»?«. <«) 



6 y V 2 dlnl s 

Note that for a scale invariant Ci s , which occurs when also the short modes are out of the Hubble 
radius at recombination, the second term in parentheses cancels and we recover the bispectrum 
computed in [1] in the squeezed limit. 

In fig.[T]we plot the angle independent part of eq. (45), i.e. when we neglect the 6 cos 29. We 



see that the final result is quite small over the whole range of Is- /nl ~ 1 - The ^-dependent 
part is larger by a factor of 6, but it will not contaminate any primordial signal, as it averages 
to zero when summed over all relative orientations of the mo des0 

Notice that there is a partial cancellation between the effect coming from the redefinition 
of the spatial coordinates and the isotropic part of lensing, as already noticed in [Tj. This is 



4 All primordial shapes, to our knowledge, have an which is independent on the relative orientation 
of the modes in the squeezed limit. 
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easy to understand as overdense regions around the surface of last scattering will have a positive 
convergence, and thus move the short-scale 2-point function towards larger angular scales; at the 



same time these regions will have $ < and this, through eq. (30), moves the 2-point function 
in the opposite direction]^] This cancellation considerably suppresses the final result, so that 
it would be quite misleading to separate lensing from the other effects. Notice that the effect 
of lensing close to the last scattering surface is included (though dwarfed by the much larger 
ISW-lensing signal) in the studies of lensing-induced bispectrum QjJJ 1201 ED 122] ; therefore, one 
must be careful not to double-count it. 



The ^-dependent term of eq. (44), which gives zero when we average over the relative orien- 



tation of the long and short modes, is entirely due to the shear (i.e. the traceless) part of lensing. 
This term gives zero when the short-scale 2-point function has Ci s = const, i.e. for a white-noise 
spectrum, with a real-space correlation which is a delta function. This is easy to understand as 
shear just remaps the points, preserving the area, so that a white noise remains such (see for 



example |22j). Also the first term in eq. (44), which includes the isotropic part of lensing and 



equation ( 42 ) , cancels for a white spectrum although we did not find any simple explanation for 
that 13 



To end this section, let us discuss an important point of eqs. (42), (43) and (44). As we 



said many times, the effect of a super-Hubble mode on the short ones simply amounts to a 
coordinate redefinition of the background. Note that the long- wavelength modes ki that we are 



5 In pQ a wrong overall minus sign in front of the two effects was introduced, when moving from real 
to Fourier space. This does not affect the partial cancellation between the two contributions. 



Notice that eq. (43 1 holds only for a scale-invariant spectrum of the long mode. Indeed, when we 



deviate from scale-invariance, the relationship between the long-mode spectrum Ci L and the correlation 
Cj^ between temperature fluctuation and lensing potential changes. In flat sky the temperature power 
spectrum in the Sachs- Wolfe limit is given by 



C lT = 



3D 2 



dk i 



(k 2 ± + ja/D3)(3-(n.-l))/2 



A 



2^ r 



D 1 



r 



(46) 



where A is the normalization of the power spectrum, D = rjo — 77* is the distance from the last-scattering 
surface and n s — 1 is the deviation from a scale-invariant spectrum. In the same limit the correlation 
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considering are already super-horizon when the short modes are quantum mechanically generated 
during inflation and they remain so until after recombination. Thus, also the primordial 2-point 
function of the short modes is rescaled by the presence of the long-wavelength modes. This 
implies that our results, eqs. (42), (43) and (44), already take into account the single-field 
Maldacena consistency relation |25[ 126] [27] , which relates the primordial 3-point function in the 
squeezed limit to the spectral tilt. Indeed, the tilt of the primordial spectrum contributes to 
the scale dependence of the temperature 2-point function above. In other words, our results 
already include the contribution from the primordial non-Gaussianity in any single-field model. 
Indeed, we checked that this is the case also when a substantial departure from scale-invariance 
in the primordial spectrum is allowed. Multi-field models of inflation can produce additional 
non-Gaussianity in the squeezed limit (see for instance [HJ E21 [33] ) and this should be added to 
our results. 



5 Comparison with numerical calculations 

Before comparing the results of the previous section with CMBquick, we would like to stress 
that the coordinate transformation derived in section [2] can be used to reproduce the evolution 
of perturbations at second order in the squeezed limit. In particular, eqs. (16), (17), (20) and 



(22) give the same evolution as that obtained from Einstein's equations when one of the two 



modes in the second-order source terms remains out of the Hubble radius. 



As an example, we consider the radiation to matter transition where e is given by eq. (12) 



and we study the second-order evolution of a small-scale mode of the potential <I> and of the radi- 
ation energy density perturbation 5 rac j in the perfect fluid approximation. Fourier transforming 



eqs. (16) and (20) (we use G = <5 ra d/4) and using eq. (12) yields 



4 2 \k s ,n) = f( V ) 



dT£\k s ,v) dT^(k s ,n) 



dlnn dinks 



(49) 
(50) 



d In r\ dinks 

where f(n) can be read from eq. (12) and by Z"W and we denote the first and second-order 



transfer functions. For a generic perturbation X these are defined by 

4%) = T«(M)C S , 



4° fa) 



between temperature and lensing potential is given by 



C, 



Tip 



3D 2 



+ OG 



v - V* 



A 



„ dk± J rh dV D* (fci+Z|/ J D2)(3-K-i))/ 2 



-A 



V7rr( 



D 



(51) 
(52) 

l-ra. 



sr (2 - ^) i\ 



We see that the relative coefficient between the two depends on the tilt 



Tiji 



3 - n s 
'l 



(47) 



(48) 



Thus, the overall coefficient of eq. (43) changes when departing from scale-invariance on large scales, 



while eq. (42) does not. This, in particular, means that the cancellation between the two effects for a 
white-noise spectrum on short scales, C/ s = const, only occurs for a scale-invariant large-scale spectrum. 
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Figure 2: Comparison of the effect of the coordinate transformation ^ on the gravitational potential 
$ and on the radiation energy density <5 ra d = 40 to the full second-order evolution in the squeezed 
limit. We consider the radiation-to-matter transition in the fluid approximation. On the left-hand panel 



compare eq. ( 49 ) (solid line) to the numerical solution of Einstein's equation (dashed line) . On the 



right-hand panel we do the same for eq. ( |50[ ). We chose fej, = 0.01?7 , k s = 10r] CQ and cos 9 = 0.5 



In fig. [2] we compare eqs. (49) and (50) with the solution of Einstein's equations at second order 
(see, for instance, pS])- We chose ki = O.Olr/" 1 , ks = lOr/' 1 and cos# = 0.5, where 9 is the 



/eq > "O -"-"'/eq 

angle between the long and the short mode and r] eq is the time of radiation-matter equality. The 
agreement is very good and gets better as we increase the ratio ks/ki. 

Let us go back to the bispectrum. As discussed before, in the squeezed limit our results 
can be also used as a consistency check of numerical Boltzmann codes at second order. Here 
we compare the bispectrum computed in the previous section with that produced by CMBquick 

m- 

Since the lensing bispectrum is a well-known result in the literature, we have chosen to remove 
the lensing effect from the code. This is easy to do, because the code parallels our calculation 
using Newtonian gauge. In this way we do not have any ^-dependence in the bispectrum and 
we avoid numerical cancellations between the isotropic part of lensing and the space redefinition 
effect. Thus, we concentrate on reproducing the bispectrum given by eq. (42). The comparison 
is direct because CMBquick computes the bispectrum using the flat-sky approximation. 

In fig. [3] we compare the bispectrum from CMBquick with eq. ( 42 ) . Solid lines represent 
our equation while dashed lines the bispectrum computed by the code. In all plots we have 
used the best-fit cosmological parameters of WMAP5 [34J. In the top panels we have fixed one 
of the multipoles to be l\ = 20 and we have plotted the bispectra as a function of the other 
two, taken to be equal h = h = I- We expect the two curves to converge in the limit I 3> 20. 
To help the comparison at low multipoles, in the bottom panels we have taken l\ = Z /10 and 
I2 = I3 = I. Here we expect the two curves to converge for small I. The code agrees reasonably 
well with our formula, for I ^ 200 in the top panels and for I 300 in the bottom panels. 
The agreement improves when we remove neutrinos, as shown on the right-hand panels. In 
this case, in order to maintain the same amount of radiation as in the left-hand panels, we 
have increased the photon temperature today to T = 3.08 K. The reason for this improvement 
is interely numerical. Neutrinos require finer integration steps and by removing them we can 
reach a better accuracy for the same computing time. There is another important source of 
inaccuracy in the results of the code. As the second-order numerical integration is particularly 
slow, we have integrated the second-order sources along the line of sight only in a finite range 
of conformal time, i.e. between 77 = 230 Mpc to rj = 350 Mpc (the peak of the visibility function 
being at ~ 280 Mpc), neglecting part of the early ISW effect. In conclusion, the code is reliable in 
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Figure 3: Comparison of eq. (42) with the bispectrum in the fiat-sky approximation produced by CM- 
Bquick in the absence of lensing. Cosmological parameters are the best-fit ones of WMAP5 [34]. Solid 
lines represent our equation while dashed lines the bispectrum computed by the code. In the top panels 
li = 20 and h = h = h m the bottom panels l\ = 1/10 and I2 = h = I- For the top panels the squeezed 
limit corresponds to / £ 300 and for the bottom panels to I ,$ 300. On the right-hand panels we have 
removed neutrinos and set the photon temperature today to T — 3.08 K to have the same radiation 
density. 



the regime that we studied and we are confident that the residual discrepancies can be removed 
by a proper optimization. Note that, in comparing the code with our analytical approach, we 
corrected some numerical inaccuracies which may affect the results of the current version of [16 . 



6 Contamination and observability 

In the flat-sky approximation, the estimator of a local primordial signal is given by |35] [36] 

If d 2 h d 2 l 2 d 2 l 3 2 r> - - B\ oc (li,l 2 ,h) /j* \ iT\ n \ / M \ 

£loc = Jf) (2^F (2^F W? { ] ( 1 + 2 + ;3) lw7 a(ll)fl(l2)fl(i3) ' (53) 

where B\ oc {li,l2,l 3 ) is the bispectrum for local non-Gaussianities when = 1 and the nor- 
malization factor that makes the estimator unbiased {{£\ oc ) = /nl) ^ s 

1 f d% d%[B loc (h J2,h)] 2 , , 

n J (2tt)2 QC h C h C h ■ {i) ' 
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Thus, if the local estimator above is applied to an arbitrary signal Bxifa, fa, fa), this will con 

of /nl g iven b y 

d 2 fad% B loc (h,fa,fa)B x (h,fa,fa) 



taminate the measurement by a value of given by 



i-loc 
/NL 



1 

M 



7T 



(2*) 



(55) 



Let us apply eq. (55) to our bispectrum (44) to compute its contamination on a measurement 
of a local signal. We evaluate the integral^] in the squeezed limit and take into account the 
triangular inequalities by choosing 2 < fa < 100, 20 fa < fa < Z max and fa — fa < fa < min(/i + 
fa,lmax)- We find 

/I- = _o. 39 f (57) 

for Z max = 2000. This value is visually confirmed by fig. [TJ where, for I ~ 2000, the plotted 
equivalent oscillates around ~ —0.4. Note that the anisotropic part of lensing, proportional 
to cos 28, does not contribute to this value because the integration over the angle averages to zero. 
However, the cancellation due to the isotropic part of lensing is crucial to obtain such a small 
result]^] As the contamination is given by the ratio between two integrals, this result is quite 
solid and changes little when we vary the range of integration keeping /max 

fixed. We conclude 

that the bias introduced by this effect on a Planck search for a primordial non-Gaussian signal 
will be negligible. The contamination increases if we take higher values of l m3uX , as confirmed 
by the behavior of in fig. [I] For instance, = —0.48 for a futuristic experiment with 

^max = 3000. 




0.001 



150 200 



500 700 

^max 



1500 2000 3000 



Figure 4: Signal-to-noise ratio for the bispectrum in eq. (44), evaluated summing over squeezed config- 

in- 7 / 2 



urations only, as a function of ? max - This roughly oscillates around 
dashed line shows S/N for a primordial local signal with = 1- 



For comparison, the 



Let us study now the observability of the bispectrum in eq. ( 44 ) . The possibility of measuring 
a bispectrum depends on its signal-to-noise ratio S/N. In the flat-sky approximation this is given 



7 In order to solve the integrals in eq. (55) it is convenient to use that 

lidh l 2 dl 2 fadfa 



d% d%F{hM,fa 



4-7T 



Hi nil + nil + mi 2 3 ) i/2 



F{fa,h,fa) 



(56) 



where F is a function of the moduli of the multipoles. 

8 For ? max = 2000, the contamination from the bispectrum in eq. (42) is /^l = 0.94, while that from 
lensing alone, eq. (43), is fj^£ = —1.33. 
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by [36] 

S\ 2 1 f d 2 l 2 d 2 h[B(h,l 2 ,l 3 )f 



NJ nj (2tt) 2 6C h C h C h 



(58) 



Plugging eq. (44) in this expression and considering all cyclic permutations we obtain the signal- 
to-noise shown in fig. [4] (solid line) as a function of the maximum multipole l max - For comparison, 
we have also plotted the signal-to-noise for a primordial local signal with /^l = 1 (dashed line). 
The signal-to-noise from eq. p4"| ) grows approximately as Z^ ax and not as Z max as one would 
naively expect by simply counting the number of available modes. Indeed, the effect of the 
exponential damping at high-/ contributes to the logarithmic derivative in eq. ( 45 ) with a term 
which increases roughly linearly in Z max - Note that, in the estimate of the signal-to-noise, the 
anisotropic part of lensing does not cancel and gives a large contribution to the signal. 

From fig. [4] we see that this effect is not measurable by Planck. Note however that we have 
integrated only over configurations in the particular squeezed limit in which one of the modes 
is out of the Hubble radius, i.e. for II < 100. Thus, this S/N represents only a fraction of what 
one would get by summing over all configurations of the full bispectrum generated by a complete 
numerical Boltzmann code. The fact that the contamination from second-order perturbations 
will not be relevant for Planck is confirmed by other studies [111 [T4"l I15j . which focused on 
particular second-order effects. 

Note that even though eq. (44) includes the primordial non-Gaussianity in the squeezed limit 
for single-field models, for the tiny value of the spectral index that we consider (n s — 1 = 0.037), 
this primordial contribution is negligible so that the result can be interpreted as being solely 
due to nonlinearities at recombination. 



7 Conclusions 

Although the calculation of the CMB anisotropics at second order is remarkably complicated, 
the CMB bispectrum in the squeezed limit, with one of the modes outside the Hubble radius 



at recombination, is given by the simple formula (44). This expression takes into account all 
the effects at recombination and must be supplemented by the additional contributions to the 
bispectrum at low redshift. As we stressed, the formula already takes into account the primordial 
non-Gaussianity produced by any single-field model of inflation, at any order in slow-roll. From 



this point of view, eq. (44) can be seen as an extension of Maldacena's consistency relation 



25"! [26| [27] in terms of observable quantities: a deviation of the temperature bispectrum from 



our eq. (44) (once low-redshift effects are properly subtracted) would rule out all single-field 
models. 

Of course, to fully exploit the forthcoming Planck data, one wants to consider also config- 



urations of the bispectrum for which (44) does not apply. In this case one has to resort to a 
numerical Boltzmann code. Our simple formula can be used as a non-trivial check of the nu- 
merics, as all the second-order effects must exactly conspire to match our simple result in the 
squeezed limit. The current version of CMBquick is able to correctly reproduce our analytic 
results, though some residual disagreement indicates that some effort should be put to improve 
the numerical reliability of the code. 

A further handle to constrain primordial non-Gaussianity will come from polarization mea- 
surements |37j . Our analysis can be extended to bispectra including polarization, along the lines 



indicated by [lj. We leave this for future work. 
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A CMB anisotropics at second order 

In this appendix we provide a second-order expression for the CMB temperature anisotropies in 
the perfect fluid approximation and assuming instantaneous recombination. We are interested 
in computing the CMB temperature fluctuations 

n, f-\ — -^obs(^) — (T bs) /_„n 
©obs(n) = -=— , (59) 

\Jobs/ 

where T \> s {n) is the observed photon temperature in the angular direction h (h 2 = 1) and (T b s ) 
is its average over the sky. We concentrate on the limit of instantaneous recombination so that 
the decoupling takes place at ci given pliysicEil time T/rec* 

For a black-body spectrum the observed 
temperature T b s (ra) is related to the one of emission T rec by Liouville's theorem; as phase space 
density is conserved in the propagation of photons, the phase space density observed in a given 
direction h is the same as at emission but with a rescaled temperature, 

T bs(n) = ^^r re c(r/rec, x vec ) , (60) 

-C'rec 

where E vec and E Q ^ S are the energies at emission and observation of a given photon. This 
equation is exact and therefore holds at any order in perturbation theory. 

We consider a flat FRW metric perturbed at second-order in Poisson gauge [38] , 

ds 2 = a 2 (r/) [-e 2 % 2 + 2uj i dx i dr ] + (e~ M ' 5^ + 7ij ) dx l dx j ] , (61) 

where w$ is transverse, d l LOi = 0, and 7^ is transverse and traceless, d l ^ij = = 7]. We only 
consider scalar primordial perturbations. In this case the vector and tensor components uii and 
jij, are only second-order quantities. 

The energy of a photon as measured by an observer with four- velocity u^ 1 , with normalization 
condition u^u^ = —1, is given by 

E = -PpvP , (62) 

where P M = dx^/dX is the four-momentum of the photon, P^P^ = 0. Thus, from eq. ( |62[ ) we 
obtain ^ 

E T ec Oobs-Po^rec) \ Pq J rec 

where we have used u° = y/l + v 2 e _ */a from the normalization condition of the four- velocity 
and the definition of v % given above, v l = ae~^u l . We have also set to zero the velocity of the 
observer, which is simply responsible for a dipole. 

The photon geodesic equation can be conveniently rewritten as [39] 



In order to compute Po we need to solve the time component of eq. (64) up to second order. 
After few manipulations we can rewrite this, up to second order, as 

1 flP " + dft - li'ijhW , (65) 



P drj 1 2 
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where we have used that h % = P l /P° = constant along the unperturbed geodesic. Upon inte- 
gration, this equation yields 



Poiv, 



'obs) 



exp [I(r] mc ,r] ohs )] , I(rj 

rec j ^obs J — 



»?obs 
'Hrcc 



dr] ($' + *'- oj'ih 1 



1 / - i ~ j 



(66) 



Note that Pn is conserved at zeroth order and we choose Pn = 1. 



As v % in eq. (63) is a first-order quantity, we need to find Pi{r] rec ) / Po(r] rec ) up to first order 

(67) 



only. Thus, we solve the spatial component of eq. (64). At first order this reads 

1 dR 



P dr] 

which can be integrated using the zeroth-order value for Pq giving 



fVoba 

PiiVrec) = PiiVobs) + dr] + 

J f?rcc 



Furthermore, eq. (66) at first order reads 

Po(r]rec) = -Po(%b s ) 



'/obs 



dr] + #') . 



(69) 



The observer photon direction is fi % = P^ hs /P^ 



obs - -Pi(Vobs)/Po(Vohs)- Here we have set to zero 
the metric perturbations at the observer position: since metric perturbations at the observer do 
not depend on the direction of observation, they can be reabsorbed in the redefinition of the 
average temperature. The two equations above can be combined to give 



Pi(Vr. 



P (Vr 



-h\l - $ rec - I' rec) + <5n*(r/ rec , r] obs ) , 5n l (r] 

rec j ?7obs ) 



'/obs 



efyVj^ + tf) , (70) 



where we have rewritten the spatial derivative in (68) as a spatial gradient orthogonal to the 
line of sight, = (Sij — hihj)dj, while the derivative along the line of sight gives a boundary 
term and an integral of <£' + 4 f/ which cancels with the last term of eq. ( 69 ) . 



Using eqs. (66) and (70), eq. (63) becomes, up to second order, 

-Ebbs 



O-obs 



exp [<J> rcc + 



I (Vrec, Vobs)] (V 1 + v mc ~ » ' ^rec + 5n(r] rec j ?]obs ) ' V r 



(71) 



Finally, plugging this expression into the expression for the temperature, eq. (60), using the 
definition of temperature perturbation, eq. (59), and expanding up to second order, we obtain 

1 1 

©obs(ra) = 6 rcc + $rec + n ■ V rec + I(r] Iec , 7/obs) + $rec@rec + T^rec ~ ^rec 

+ (fl ■ V TCC + 7(??rec, ??obs))(©rec + ^rec + fl ■ V rec ) - 5n(r] rec , T] ohs ) ■ V rcc + -I(j] rec ; 'Hobs ) 

(72) 

It is straightforward to verify that this expression agrees with the one give in [3]. 

For the sake of clarity we choose to simplify this expression even further by assuming pure 
matter dominance and absence of anisotropic stresses. In matter dominance I(i] vec , i] obs ) is only 
second order since at first order <E> = = const. Moreover, this second-order contribution (which 
includes the Rees-Sciama effect, vector and tensor modes) does not affect the CMB bispectrum 
in the squeezed limit [I], so that we can neglect all terms involving / in eq. (72). Furthermore, 
as v vanishes in the large scale limit, we can neglect in eq. (72) all terms which are quadratic in 
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v. It is easy to realize that also the term 
be neglected^ 



-5n ■ v Iec vanishes in the squeezed limit and can thus 
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This leaves us with 



©obs(^) =[Q + $ + h-V + <S>0 + $ 2 /2 + (9 + $)n • v] (?7rec,^rec) • 



(74) 



Note that the right-hand side of this equation depends on the physical time r/ rec and position 
x Tec of photon emission, which are perturbed quantities. Expanding around the unperturbed 
emission time and position gives rise to other second-order effects, well described in [40J. We 
discuss these effects in more details in the main text. 



http://www.rssd. esa.int/index.php?project=Planck 
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